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The scenario of SO (10) grand unified theory (GUT) proposed by one of the authors is 
extended to Ee unification. This gives realistic quark and lepton mass matrices. In the 
O . neutrino sector, the model reproduces the large mixing angle (LMA) MSW solution as well 

' as the large mixing angle for the atmospheric neutrino anomaly. In this model, the right- 

handed down quark and the left-handed lepton of the first and second generations belong 
to a single multiplct 27. This causes natural suppression of the flavour changing neutral 
current (FCNC). 

§1. Introduction 

00 ' 

There are strong reasons to believe in the validity of grand unified theories 
(GUT) EP , in which the quarks and leptons are unified in several multiplets in a simple 
gauge group. They explain various matters that cannot be understood within the 
standard model: the 'miracle' of anomaly cancellation between quarks and leptons, 
the hierarchy of gauge couplings, charge quantization, etc. The three gauge groups 
in the standard model are unified into a simple gauge group at a GUT scale, which 
is considered to be just below the Planck scale. On the other hand, the GUT scale 



destabilizes the weak scale. One of the most promising ways to avoid this problem 
is to introduce supersymmetry (SUSY). However, it is not easy to obtain a realistic 
SUSY GUT.U First, it is difficult to obtain realistic fermion mass matrices. Also 
unification of quarks and leptons puts strong constraints on the Yukawa couplings. 
Finally, one of the most difficult obstacles is the "doublet-triplet (DT) splitting 
problem." 

There have been several attempts to solve the DT splitting problem, u^o 1 Among 
them, the Dimopoulos-Wilczek mechanism is a promising way to realize DT splitting 
in the 50(10) SUSY GUT.ft'B'B-EP 

Concerning the fermion masses, recent progress in neutrino experimentstP pro- 
vides i mpo rtant information on family structure. There are several impressing 
works Q^'EJ'tiP' in which the large neutrino mixing angles are realized within 

GUT framework. It is now natural to examine SO(W) and higher gauge groups 
because they allow for every quark and lepton, including the right-handed neutrino, 
to be unified in a single multiplet, whic h is important in addressing neutrino masses. 

Recently, one of the authors (N. M.)c3* proposed a scenario of SO(W) grand uni- 
fied theory (GUT) with anomalous U(1)a gauge symmetry, which has the following 
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interesting features: 

1. The doublet-triplet (DT) splitting is realized using the Dimopoulos-Wilczek 
mechanism. 

2. Proton decay via the dimension-five operator is suppressed. 

3. Realistic quark and lepton mass matrices can be obtained in a simple way. In 
particular, in the neutrino sector, bi-large neutrino mixing is realized. 

4. The symmetry breaking scales are determined by the anomalous U (1)a charges. 

5. The mass spectrum of the super heavy particles is fixed by the anomalous U(1)a 
charges. 

As a consequence of the above features, the fact that the GUT scale is smaller 
than the Planck scale leads to modification of the undesirable GUT relation between 
the Yukawa couplings = y s (and also y e = y^) while preserving the relation 
Ut = Ub- Moreover, it is remarkable that the interaction is generic; all the interactions 
that are allowed by the symmetry are taken into account. Therefore, once we fix 
the field contents with their quantum numbers, all the interactions are determined, 
except coefficients of order 1. 

The anomalous U(1)a ga uge symmetry,© whose anomaly is cancelled by the 
Green-Schwarz mechanism, t3 plays an essential role in explaining the DT splitting 
•hanism at the unification scale as well as in reproducing Yukawa hierarchies. 
I8|),|i9p Also, bi-large neutrino mixing is naturally obtained by choosing the 10 
representation with an appropriate U(1)a charge, in addition to the three family 16 
representations. This anomalous U(1)a is a powerful tool not only to reproduce DT 
splitting but also to determine the GUT breaking scales. 

This paper aims to show further that the above 50(10) model is naturally 
extended to Eq GUT, in which the additional field 10 of SO(W) is included in a 
chiral multiplet 27 of Eq. In order to realize this scenario, it is important to introduce 
the concept of "twisting family structure" in the Eq unified model. c3 

Under the 50(10) group, we know that 10 and 5 of SU(5) are combined into 
16 of 5*0(10). Usually, each family belongs to 16. In this framework, however, it is 
not easy to reproduce the large Maki-Nakagawa-Sakata (MNS)ESP mixing and small 
Cabbibo-Kobayashi-Maskawa (CKM)@> mixing. A promising way to reproduce this 
is to introduce other multiplets, 10 of 5O(10), in addition to the usual 3 x 16 
multiplets.0* Since 10 of SO (10) is decomposed into 5(10) and 5(10) of SU(5), 
one of the fields 5(16) can be replaced by this 5(10). Such a replacement is essential 
to reproduce large MNS mixing, preserving small CKM mixing. In the case of Eq, 
16 and 10 of 5O(10) are naturally included in a single multiplet 27 of Eq. The Eq 
model automatically prepares such a replacement, as we see in the next section. We 
call the mechanism responsible for this the "twisting mechanism." This gives us a 
strong motivation to examine Eq unification. 

It is interesting that the above desired scenario in 5O(10) unification can be 
extended to Eq unification while keeping the desirable features of 5O(10) unification. 
In this paper we focus on the extension of the matter sector to Eq unification, leaving 
discussion of DT splitting to a separate paper.© The extension of DT splitting to Eq 
unification is non-trivial. Moreover, we show that the condition for the suppression 
of the flavour changing neutral current (FCNC) is automatically satisfied. This is 
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essentially caused by the twisting mechanism and the unification of the matter fields 
into a single multiplet 27, which guarantees that 5(16) and 5(10) have the same 
anomalous U(1)a charge. Then it can happen that the charge of the first generation 
of 5 becomes equivalent to that of the second generation of 5. This weakens the 
severe constraint resulting from the FCNC. It is interesting that the selection of 
the anomalous U(1)a charge to realize bi-large neutrino mixing angles automatically 
causes the above mensioned FCNC suppression. 

In section 2, we briefly review the twisting mechanism and classify the patterns 
of the massless modes of the 5 fields. In section 3 and the Appendix, we discuss 
how the SUSY vacua are determined in the anomalous U(1)a framework. In section 
4 and 5, we study the realistic quark and lepton mass matrices in Eq unification, 
where in the neutrino sector, bi-large neutrino mixing angles are naturally realized. 
In section 6, we examine the effect of SUSY breaking. Specifically, we are able to 
automatically obtain a condition for suppression of flavour changing neutral current 
process (K°K° mixing) in Eq unification. 

§2. Twisting in Eq unification 

Let us first recall the twisting mechanism, which has been proposed by one of the 
authors (M. F3. ) . HiP The twisting family structure arising through this mechanism is 
peculiar to the Eq unification model, and here we explain how it arises. In the case 
of Eq, 16 and 10 of SO(10) are naturally included in a single multiplet 27 of Eq. 
The fundamental representation of Eq contains 16 and 10 of 50(10) automatically: 
Under Eq d 50(10) D SU(5), we have 

27 -> [(16, 10) + (16, 5) + (16, 1)] + [(10, 5) + (10, 5)] + [(1, 1)], (2-1) 
16 10 1 

where the representations of 50(10) and SU(5) are explicitly denoted. As we have 
already seen, the Eq model naturally possesses the freedom to replace matter fields 
(16, 5) by (10, 5). So here let us explain how the twisting family structure arises in 
the Er unification. In order to do this, it is enough to introduce the following Higgs 
fieldsjil which are necessary to determine the mass matrices of matter multiplets 
#i(27), whose U(1)a charges are denoted as ipf^ii = 1,2,3): 

1. A Higgs field that breaks Eq into 5O(10): <2>(27) (($(1, 1)) = v). 

2. A Higgs field that breaks 5O(10) into SU(5): 0(27) ( (0(16, 1)) = v'). 

3. A Higgs field that includes the Higgs doublets: H(27). 

Throughout this paper we denote all the superfields with uppercase letters and their 
anomalous U(1)a charges with the corresponding lowercase letters. Assigning nega- 
tive R-parity to the ordinary matter ^(27), as usual, and using a field with charge 

*' Note that the additional Higgs fields #(27) and (7(27) are required to satisfy the Z)-flatness 
condition of E% gauge theory, and an adjoint field ,4.(78) is required to break the GUT gauge group 
into the standard gauge group. In order to realize doublet-triplet splitting, the actual breaking 
pattern-must be E 6 -> SO (10) -» SU(3) C x SU(2) L x SU(2) R x U(1)b-l -» SU(3)c x SU(2) L x 
U(1)yMM> 

**•' We assume that > ip2 > i>3 
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— 1, the U(1)a invariant superpotential for low energy Yukawa terms becomes 

Wy = (— ) ##;il, (2-2) 

where we omit coefficients of order 1, and for the above we assume that ipi+ipj+h > 
for each i,j pair so that there appears no SUSY zero.0 After obtaining non-zero 
VEV (0) = XMp (A ~ 0.2) through the D-flatness condition of the anomalous 
U(1)a gauge symmetry, a hierarchical structure of the Yukawa couplings is realized. 
Since we need 3 x [10 + 5 + 1] in SU (5) representations for three families, among 
the above three 27 fields, three pairs of (5, 5) must become heavy.El] Indeed, the 
Higgs fields <P and C can yield such masses: The superpotentials, which give large 
masses for (5, 5) pairs, are 

w={ — ) ##,<? + [ — ) ##,#, (2-3) 



which we analyse here to see how 5 fields acquire large masses. The VEV (#(1, 1)) 
v gives the 3x3 mass matrix of #i(10, 5)3^(10, 5) pairs, 

#i(10,5) # 2 (10,5) # 3 (10,5) 

!Pi(10,5)/ A 2 * A*+^ 2 A^ 1+ ^ 3 \ 
{Mij)= #2(10,5) A^ 1+ ^ 2 A 2 ^ 2 A^ 2 +^ AV (2-4) 
3>3(10,B) \ \^+^ A^ 2+ ^ 3 A 2 ^ 3 J 

while the VEV (C(16, 1)) = v' gives the mass terms of #(16, 5) and #,(10, 5), 

#i(16,5) # 2 (16,5) # 3 (16,5) 



#i(10,5) 

(mm= #2(10,5) 



/ ^2l/>l ^l/>l+</>2 \1pl+f3 

\lpl+lp2 \2lP2 _\</>2+V>3 | \Cy\ (2-5) 

\tpl+lp3 \4>2+4>3 \2i>3 



#3(10,5 

Then, the full mass matrix is 

#i(16,5) # 2 (16,5) # 3 (16,5) #i(10,5) #2(10,5) #3(10,5) 

#i(10,5) / A 2 ^ 1+r yPi+^+r \4>i+4>3+r ^2-01 \4>i+4>2 \i>i+-4>3 \ 

#2(10,5) X^l+i'i+r )2ip 2 +r yl>2+i>3+r yl>l+i>2 \2ip2 yl>2+i>3 \\ ( t> l 

#3(10, 5) I A^ 1+ ^ 3+r A^ 2+ ^ 3+r A 2 ^ 3+r XPi+^3 \Tp2+1p3 \2lp3 J 

(2-6) 

where we have defined the parameter r as 



AV 



*' Note that if the total charge of an operator is negative, the U(1)a invariance forbids the 
existence of operators in the action, since the field with negative charge cannot compensate for 
the negative total charge of the operator (the SUSY zero mechanism). 

**' The possible right-handed neutrino modes ^(16,1) and 1) also acquire large masses, 

but here we concentrate on the family structure of 5. 
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which we use frequently in the following discussion. Note that some of the matrix 
elements become zero if the index becomes negative (the SUSY zero). For the mo- 
ment we assume that no such zero appears in the superpotential. In general, it is 
seen that we have three massless modes out of the six 5 fields by solving the above 
3x6 matrix. However, since the matrix has hierarchical structure, we can easily 
classify the cases. 

1. Under the condition that we have no SUSY zeros, it is evident that the largest 
mass is either M33 or M33, whose ratio is M33/M33 = A r . 

2. < r case: In this case M33 is larger than M33, and the pair ($3(10, 5), $3(10, 5)) 
is heavy. Next, compare M 23 and M 22 , whose ratio is M 23 /M 22 = \ r+ ^-^ 2 . 
Thus, there are several cases, depending on r and ^3 — V'2- 

3. < r < ifj 2 - ip3- In this case, M 23 > M 22 , so that the pair $3(16, 5)# 2 (10, 5) 
becomes heavy, and at the same time the pair ^2 (10, 5)#i(10, 5) obtains a large 
mass, because M' l2 /M 12 = X r < 1. Oq(10,5), #i(16,5) and #2(16,5) are left 
massless. This case is denoted (1,1', 2). [In this paper, massless mode whose 
dominant component is #j(16,5) (^(10,5)) is simply denoted by 

4. (0 <)ifj 2 - Y>3 < r: The pair # 2 (10, 5)# 2 (10, 5) becomes heavy. Further, this 
case is divided into two cases, according to the sign of r + ^3 — ipi ■ 

5. ^2-^3 <r < Y>i-Y> 3 : In this case, M' l3 is larger than M u . Thus # 3 (16, 5)#i(10, 5) 
becomes heavy, and this case also becomes the case (1, 1', 2). 

6. (^2-^3 <)ipi-4>3 < r: In this extreme case, M' n < M n , and thus #1 (10, 5)#i(10, 5) 
becomes heavy. Hence all the <Pj(10, 5) are heavy states and the ^(16, 5) are 
massless modes. This corresponds to the situation in which the three massless 

5 fields (quarks and leptons) belong to ^(16,5). This is just the case usu- 
ally adopted in the 5O(10) model. We call this case that of "parallel family 
structure." We denote this case simply as (1,2,3). 

7. r < case: This case is easily classified just replacing the 10 representation 
with the 16 representation. 



Thus we can classify all the cases as follows: 

1. Vi-V'3<r: (1,2,3) type. 

2. 0<r<^i-^ 3 : (1,1', 2) type. 

3. Vs-V'i < r < : (1,1', 2') type. 

4. r <^ 3 -V>i : (l',2',3') type. 

If we use SUSY zero coefficients, various types of massless modes can be realized. 
For example, if V'l + ^3 + 4> < 0, SUSY zeros appear, and the Yukawa terms #3^^ 
(i = 1, 2, 3) are forbidden. Hence the mass matrix M becomes 



and the massless mode #3(10,5) does not mix through non-diagonal mass matrix 
elements with any other 5 field. We call such a massless field an "isolated" field. 
There are various different patterns of massless modes containing the "isolated" 




(2-8) 
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fields. For example, if the conditions 2^2 + </> > 0, 2^3 + c > and X 2 ^ 1+ ^v > 
\ipi+if>2+c v ' are sa ti s fied in addition to the above condition + ^3 + 4> < 0, we have 
the pattern (1,2,3'), i.e., 

/ 5i \ / #i(16,5) + .. • \ 

5 2 = # 2 (16,5) + ••• , (2-9) 
\ 5 3 J \ #3(10,5) J 

which has been adopted in Ref. [ll]). Note that 53 does not mix with any other 
states (an isolated field). 

In addition to the mixing of the matter content, the massless Higgs doublets 
itself can in principle be mixed as 

H(S) = #(10,5)cos6> + F(16,5)sin<9, (2-10) 

which is also determined by obtaining a whole mass matrix of the doublet Higgs fields. 
Note that the Yukawa couplings #(16, 10)#,(16, 5)H(5) (#(16, 10)# i (10, 5)iT(5)) 
are proportional to cos 9 (sin 9). 

§3. Features of the vacua in the U(1)a framework 

In this section, we explain how the vacua of the Higgs fields are determined by 
the anomalous U(1)a quantum numbers. Hj) 

First, the VEV of a gauge invariant operator with positive anomalous U(1)a 
charge must vanish. Otherwise, the mechanism of the SUSY zero does not work, since 
such a VEV can compensate for the negative U{1)a charge of the term. Generically, 
such an undesired vacuum is allowed, but as is shown in the Appendix, in such a 
vacuum, the Froggatt-Nielsen mechanism c3> does not operate. Therefore we are not 
interested in such a vacuum. Here, we simply assume that we are in the vacuum 
where the SUSY zero and Froggatt-Nielsen mechanism operate, namely any VEV of 
a gauge invariant operator with positive anomalous U(1)a charge vanishes. 

Next, we show that the VEV of a gauge invariant operator O is determined by 
its U(1)a charge o as (O) = \~° if the i^-flatness condition determines the VEV. 
For simplicity, we examine this relation using singlet fields Zi with anomalous U(1)a 
charge Z{. The general superpotential is written 

W = Y, x * lZ i + lZ ^ Zl+Zj ZiZj + ■■■ (3-1) 

= E^ + E^ + ---> ( 3 - 2 ) 

i i,j 

where Zi = X Zi Zi. The equations for the F-flatness of the Z{ fields require 

A ^(l + £z j + ...) = 0, (3-3) 

j 

which generically lead to solutions Zj ~ 0(1) so that (Zi) ~ A -2 ', as stated above. 
Note that the Froggatt-Nielsen structure of Yukawa couplings, A^ i+ ^ +ft- #j# ? lf, is not 
changed by the interactions W = A^ + ^ +h+Zfc Z fc ## 7 i7 with the VEVs (Z k ) ~ A~ 2fc . 
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If an adjoint field A possesses a non-zero VEV by the F-flatness condition, this 
VEV is determined as (A) ~ A~ a , because A 2 can be gauge invariant. Suppose 
that, in addition to and C, there are 3(27) and (7(27). Since <P<P is also gauge 
invariant, the VEV of the operator is given by ~ \-($+<l>) if it is determined by 

the F-flatness condition. The D-flatness condition of the E% gauge group requires 

|(<p)| = |(^)| ^\-(4>+<t>)/i. (3.4) 

Note that these VEVs are also determined by the anomalous U(1)a charges, but they 
are different from the naive expectation (<P) ~ X~^. This is because the D-flatness 
condition plays an important role in fixing the VEVs. The VEVs of (7 and are 
also determined by the anomalous U(1)a charges as 

| = | ~ A~ (c+C " )/2 . (3-5) 

By the above argument, it is found that v and v' are determined by the anomalous 
U(1)a charges. Therefore, the massless modes of 5, which are determined by the 
twisting mechanism, are also determined by the anomalous U(1)a charges. 

§4. Quark and lepton masses in Eq unification 

Now let us consider a simple model in which realistic mass matrices of quark 
and lepton are obtained. Consider the following minimal matter content and Higgs 
chiral fields. Here, in addition to R-parity, we introduce Z2 parity, which plays 
an important role in solving the DT splitting problem, as explained in separate 
papers.®'© 

1. Matter multiplet (odd R-parity): ^(27,+) i = 1,2,3. 

2. Higgs field which breaks E 6 into 50(10): £(27,+), <£(27,+), (<f>) (=(#)). 

3. Adjoint Higgs field A(78,-), whose 50(10) component A(45) breaks 50(10) 
into SU(3) C x SU(2) L x SU{2) R x U{l) B -L by the VEV (A(45)) B _ L = ir 2 x 
diag(V, V, V, 0, 0). (This Dimopoulos-Wilczek form of the VEV plays an impor- 
tant role in solving the DT splitting problem.) 

4. Higgs field which breaks SU(2) R x U(1) B -l into U(1) Y - C(27,+), 0(27,+) 
by developing (O) (=(0». 

5. Higgs field which contains usual SU(2)l doublet: H(27, +). 

In the above, the signature db indicates the Z2 parity of the fields. Here we have 
introduced Higgs fields H in addition to the other Higgs fields <P, 3, C and O, but it 
might be the case that the Higgs doublet can be a part of a component of the other 
Higgs fields <P and/or O. Even in that case, the following argument can be applied 
by taking h = (j) or h = c. 

In the following, we take the U(1)a charges of the matter fields ^ as ipi = 
3 + n, ip2 = 2 + n and ^3 = n i which have been determined in previous papers to 
be consistent with the up-type quark masses and mixings. Then the top Yukawa 
coupling of order 1 determines the anomalous ^(l)^ charge of the Higgs field H as 
h = -2n. 

Also, in this paper, we assume that the mixing angle sin 6 [defined in Eq.flTljj)] 
is zero, i.e., the down-type Higgs is purely H(10,5). This assumption makes the 
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following analysis much simpler. Of course, once we determine the model that realizes 
the DT splitting, the Higgs mixing angle is also determined by the anomalous U(1)a 
charges. We shall discuss this point in a separate paper. c3 Actually, we find various 
DT splitting models that give sin# = 0. 

Now the Yukawa couplings are obtained by Froggatt-Nielsen mechanism! 2 ^ as 

W = X^ + ^ +h ^jH, (4-1) 

where the mass matrix of the up quark sector is uniquely determined, since we have 
already fixed the £7(1) charges of the fields #i(27): 

#i(16,10) # 2 (16,10) #3(16,10) 

#i(16,10) / A 6 A 5 A 3 \ 

M„= #2(16,10) A 5 A 4 A 2 (#(10,5)). (4-2) 

#3(16,10)^ A 3 A 2 1 J 

The twisting mechanism discussed in the section 2 causes the down quark mass 
matrix to differ from that of the up quark. 

We examine the massless modes of 5 in the following under the assumption 
sin 9 = 0. Note that in such a situation, component fields #j(10, 5), have no Yukawa 
couplings, because the Yukawa terms #i(10, 5)#(16, 10)il(10, 5) are forbidden by 
SO (10) gauge symmetry. This excludes the cases that include isolated #j(10,5) 
fields. Moreover, the cases that include isolated #i(16, 5) fields should be discarded, 
since we cannot obtain large neutrino mixing angles in such cases. With the charge 
assignment ip\ = n + 3, ip2 = n + 2 and -03 = n, the classification discussed in 
section 2 applies even to cases in which SUSY zeros appear, provided that there are 
no isolated fields. Among the options (1, 2, 3), (1, 1', 2), (1, 1', 2') and (1', 2', 3'), only 
the option (1, 1', 2) gives realistic quark and lepton mass matrices. 

Let us examine an case of (1,1', 2), i.e., < r < 3, with the constraints that 
forbid the existence of an isolated state, < ipi + V>3 + c and < + ip3 + <fi. Here, 
the parameter r, which has been defined by A r = X c v' / (\^v), is given by 

r = l[c-c -(«/>-£)], (4-3) 

because the VEVs v and v' are fixed by the anomalous U(l) a charges. Note that 
even if we take the anomalous U(1)a charges as integers, r can be a half-integer. 
This fact plays an important role in realizing the bi-large neutrino mixing angles, 
as we see in the next section. With this case, we investigate which type of mixing 
pattern of 5 fields can reproduce the bi-large neutrino mixing. In order to see this, 
let us consider the cases (5i,52,5s) = (1,1', 2) and (5i,52,53) = (1,2,1') as phe- 
nomenologically viable patterns of the massless three fields (5i,52, 53). Note that 
the correct expression of the massless states 5j at low energy are obtained as mixed 
states of #j by solving the mass matrix of Eq.( |2-6p , It should be remarked that 
#i(10, 5) itself does not have a Yukawa coupling, and therefore the field 1' really can 
have a Yukawa coupling only through the mixing with #(16, 5). In order to obtain 
the exact mass matrix for down quarks as well as leptons, we should take account of 
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the mixing effects from the non-dominant states. We first fix the three bases of the 
massless modes (5i,52,53) to (<Zq(16, 5), 3*1(10, 5), ^(16, 5)). On this basis, we can 
estimate the order of mixing parameters with the heavy states ^(16, 5), ^(lO, 5) 
and <F 3 (10, 5) as 

5t = £q(16, 5) + \^-^& 3 {16, 5) + A^-^ +r ^(10, 5) + A^-^+^lO,^) 

5 2 = ^i(10, 5) + X^-^-^W,^) + A^-^^ 2 (10, 5) + A^ 1 ^ 3 <F 3 (10 J 5)(4-5) 

5 3 = ^(16, 5) + A^ 3 <? 3 (16,5) + A r ^ 2 (10, 5) + A^ 3+r <F 3 (10, 5), (4-6) 

where the first terms on the right-hand sides are the main components of these 
massless modes, and the other terms are mixing terms with heavy states, \P 3 (16, 5), 
^2(10,5) and ^3(10,5). The order of these mixing parameters can be estimated 
by the ratios of the relevant mass matrix elements. For example, the ratio of the 
mass matrix element M' kl = A^ 1+ ^ + V to M' k3 = A^ 3+ ^ + V becomes M' kl /M' k3 = 
\tpi-ip3^ w hich appears in the coefficient of the second term of Eq. (4-4). Note 
that this ratio is independent of the parameter ip k . Similarly, the ratio of M' kl = 
A^+^+V to M k3 = Xfo-Wk+iv becomes M' kl /M k3 = A*-^ 3+r , which appears in 
the coefficient of the third term in Eq. (4-4). 

The mass matrices of the down-type quark and charged lepton can be obtained 
from the above mixing pattern by introducing the RGE factor r}~ 1 ~ 2-3. We then 
have 



£"i(16,10) 

M D = M^rf 1 = 1^2(16,10) I A 5 A 5 ~ r A 4 | (#(10, 5)) , (4-7) 
#3(16,10) 

which corresponds to the case (1,1', 2), for which 3 — r > 2 — > 1> r.O Note 
that in Eq. (4-5), the main mode of 52 is 3q(10,5), which has no Yukawa coupling 
to ,ff(10,5). Therefore the contribution from the mixing term A^ 1_ ^ 3_r !p3(16, 5) 
determines the order of the Yukawa couplings. On the other hand, the main modes 
of 5i and 53 determine the order of the Yukawa couplings, while the contribution of 
the mixing terms is of the same order. 

Now that weJbave the mass matrices for up and down quarks, we can estimate 
the CKM matrixS as 





C/ckm = A 1 A 2 , (4-8) 



which is consistent with the experimental value if we take A ~ 0.2. Since the ratio 
of the Yukawa couplings of top and bottom quarks is A 2 , a small value of tan/3 = 



*' If 1 < r, the second family should be exchanged with the third family (the case (1, 2, 1')). 
**■' Strictly speaking, if the Yukawa coupling originates only from the interaction (^^), the mixing 
involving the first generation becomes too small, due to a cancellation. In order to obtain the 



expected value of the CKM matrix as in Eq. (4-5), non-renormalizable terms, for example 'Pi'PjHCC, 
must be taken into account. 
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(#(10,5)) /(#(10,5)) ~ m t /m b ■ A 2 is predicted by these mass matrices. The 
Yukawa matrix for the charged lepton sector is the same as the transpose of at 
this stage, except for an overall factor 77 induced by the renormalization group effect. 

§5. Bi-large neutrino mixing in Eq unification 

Now we treat the neutrino masses and mixing. In order to do this, we must 
estimate the mixings in the neutrino mass matrix, since the Maki-Nakagawa-Sakata 
(MNS) matrixEP is given by 

C/mns = UiUt, (5-1) 

with the unitary matrices U\ and U v that make the matrices Ue(M^Me)U e and 
U*M v Ul diagonal. The matrix M v is the Majorana mass matrix of the light (almost) 
left-handed neutrinos, which is obtained from the Dirac masses and right-handed 
Majorana masses. First, the Dirac neutrino mass matrix is given by the 3x6 matrix 

#i(l,l) # 2 (1,1) #3(1,1) #i(16,l) # 2 (16,1) # 3 (16,1) 

5! / 
5 2 
5 3 



A r+6 


A r+5 


A r+3 


A 6 


A 5 


A 3 


A 6 


A 5 


A 3 


A 6-r 


A 5-r 


A 3-r 


A r+5 


A r+4 


A r+2 


A 5 


A 4 


A 2 



or we simply express it as 

m n = ( x r+2 a 2 ) 



(#(10, 5)) 77, 
(5-2) 

/A 4 A 3 A 

A 4 ~ r A 3 " r X 1 - 1 ' I (#(10, 5)) 77. (5-3) 



V 



A 3 A 2 1 

The right-handed Majorana masses come from the interaction 

A^ + ^ +2 ^## i ## + A^ + ^+ 5+ ^#i#,#<5 + X^ + ^ +2 %WjCC. (5-4) 

Then, the 6 x 6 matrix for #;(l,l),i = 1,2,3, and # fc (16,l),fc = 1,2,3, the right- 
handed neutrinos, is expressed as 

M R = A^+* +2 ^#(l, 1)#,(1, 1) (#) 2 + A^ + ^ +5+ ^#(l, l)# fc (16, 1) (#) (C) 
+A^+^+ 2 %(l6, l)# m (16, 1) (C) 2 (5-5) 

/ A 6 A 5 A 3 



A 2n 



^-<t> A (</»-0+c-c)/2 
A (0-0+E-c)/2 A c-c 



A 5 A 4 A 2 I , (5-6) 
^ A 3 A 2 1 

from which the neutrino mass matrix is found using the seesaw mechanismB to be 

(#(10,5)>V, (5-7) 

A X 1 - 7 - 1 J 
where we have used the relation (|4-3[). 





f 


A 2 


A 2 ~ r A 


M v = M N M^ l Mjj = A 4 - 2 "+ c - c " 




A 2-r 


A 2-2r A l-r 




\ 


A 


A 1 ^ 1 
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Combining the charged lepton sector from the previous section and neutrino 
sector from above, we finally obtain the Maki-Nakagawa-Sakata matrix:cJ 



Umns 



A' 



/ 1 

A r 1 
v A A*~ r 



A 

A l-r 
1 



(5-9) 



Recent experiments on atmospheric neutrinos have suggested a very large mixing 
angle between second and third generations, and thus r = 1/2,1 may be realistic 
[for the case of (1,1', 2), i.e., r < 1].S It turns out that r = 1/2 actually leads 
to bi-largeneutrino mixing angles, which are examined within the SO (10) model in 
Ref. |l4l).l***)| Indeed if we take r = 1/2, namely, 



+ 1, 



the MNS matrix is given by 



Umns 




(5-10) 



(5-11) 



which gives bi-large mixing angles for the neutrino sector, since A 1 / 2 ~ 0.5. At the 
same time it predicts V e ^ ~ A. It will be interesting to see if future experiments find 
evidence just below the CHOOZ upper limit V e ^ < 0.15.® For the neutrino masses, 
the model predicts rn Uii /m UT ~ A, which is consistent with the experimental data: 
1.6 x KT 3 eV 2 < Aml tm < 4 x 10" 3 eV 2 and 2 x 10" 5 eV 2 < Z\m 2 olar < 1 x 10~ 4 eV 2 , 
which is the allowed region for the most probable MSW solution for the solar neutrino 
(LMA).i) 

If we enforce the condition 



2n - 10 - I, 



the neutrino mass matrix is obtained as 

/ A 2 A 2 -'' 



My = 



A 2-r A 2-2r A l-r | (^(IQ, 5 )) 2 r? 2 , 

y a x 1 ^ i 



(5-12) 



(5-13) 



*' In the case r > 1 (1, 2, 1'), we obtain 



1 
A 



A A r 
1 A r_1 



(5-8) 



**•* In the case of (1,2, 1'), the parameter value r = 3/2 may yield a prediction consistent with 
the large mixing indicated by atmospheric neutrino experiments. 

***) When r = 1, the fermion mass matrices become of the "lopsided" type. This would seem to 
give a small mixing angle solution for the solar neutrino problem. However, recently it has been 
pointed out that taking account O(l) coefficients, lopsided-type mass matrices can give even large 
mixing angle solutions for solar neutrino problem. 
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where we have used the relation ( 5- 10| ) . From the above equation, we obtain 

'(10,5)) ; 

m Vr M P 



y = r«. (544) 



We are supposing that the cutoff scale Mp is in the range 10 16 GeV < Mp < 
10 20 GeV, which allows -2 < I < 2. If we choose 1 = 0, the neutrino masses are given 
bym UT ~ A~ 5 (H(10,5)) 2 T] 2 /M P ~ m v J\ ~ m u JX 2 . If we take 77 (#(10, 5)) = 100 
GeV, M P ~ 10 18 GeV and A = 0.2, then we get m Vr ~ 3 X 10~ 2 eV, ~ 6 x 10~ 3 
eV and m Ue ~ 1 x 10~ 3 eV. From such a rough estimation, we can obtain values 
that are nearly consistent with the experimental data for atmospheric neutrinos and 
we can also obtain a large mixing angle (LMA) MSW solution for the solar neutrino 
problem. c3> This LMA solution for the solar neutrino problem gives the best to the 
present experimental data.l 2 ^ 

Finally, we would like to make a comment on an interesting feature of this 
scenario, which is also seen in the 50(10) model. HIP In addition to Eq, fl4T|) , the 
interactions 

\^+^+ 2a + h ^A%H (5-15) 

also contribute to the Yukawa couplings after A develops a non-vanishing VEV. 
Here, only A 2 appears because of its odd Z2 parity. Since (A) is proportional to the 
generator of B — L, the contribution to the lepton Yukawa coupling is nine times 
larger than that to the quark Yukawa couplings. If we set a = —2, the additional 
matrices are 

AM U V 2 I \ " '" \ AM,, 





A 2 





A 




A 





1 


^5,16) 



















(5-17) 



(#(10,5)) 4 | n n n )' (#(10,5)) 

AM e 
(#(10,5)) " 

Note that the additional terms contribute mainly to the lepton sector. It is interesting 
that this modification essentially changes the mass eigenvalues of only the first and 
second generation. Hence it is natural to expect that a realistic mass pattern can be 
obtained by this modification: It changes the unrealistic prediction = m s at the 
GUT scale while preserving the beautiful prediction nib = m T at the GUT scale (the 
GUT relation) This enhancement factor of 2 — 3 of can be enough to improve 
the unwanted situation of the lepton quark relation in the second family. 

Remarkably enough, this charge assignment of A determines the scale of (A) 
as ~ A 2 . This strong correlation of the unification scale, which is a bit smaller 
than the Planck scale, and the improvement of the undesired GUT prediction = 
m s is indeed a consequence of U(1)a- It is also interesting that the SUSY zero 



Strictly speaking that are forbidden by the SUSY zero mechanism are generically induced by 
integrating out heavy fields that are introduced to solve the DT splitting problem. These terms 
may give a small correction to the GUT relation mj = m T . 
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plays an essential role again. When z, z > —4, the terms X^ 1+ ^^ +a+z+h Z^A^jH + 
ypi+ipj+2z+h z^q/^jff a i so contribute to the fermion mass matrices, though only to 
the first generation. 

§6. SUSY breaking and FCNC 

Finally, we discuss SUSY breaking. Since we should choose the anomalous U(1)a 
charges dependent on the flavour to produce the hierarchy of Yukawa couplings, 
generically non-degenerate scalar fermion masses are induced through the anomalous 
U(1)a -D-term.0 Various experiments on FCNC processes give strong constraints 
on the off-diagonal terms A in the sfermion mass matrices due to the fact that the 
flavour-changing terms appear only in the off-diagonal parts of the sfermion prop- 
agators, as seen in Ref. ^). The sfermion propagators can be expanded in terms 
of 5 = A /m 2 , where rh is the average sfermion mass. As long as A is sufficiently 
smaller than m 2 , it is enough to take the first term of this expansion, and then, the 
experimental information concerning FCNC and CP violating phenomena is trans- 
lated into the upper bounds on these (5fj)xY, where F = U, D, N, E, the chirality 
index is X, Y = L,R, and the generation index is i,j = 1,2,3. For example, the 
experimental value of K — K° mixing gives 



\M%)ll(S&)rr\ < 2-8 x 1(T 3 (^f 1 ) , (6-1) 
|Re(5f 2 ) Li |, \Re(5°) RR \ < 4.0 x 1(T 2 (^f 1 ) , (6-2) 



with rhn the average value of the squark masses.Liy The /i, — > ey process gives 



\(Sf 2 ) RR \ < 3.8 x 10- 3 (^P) , (6-3) 

where rhi is the average mass of the scalar leptons. In the usual anomalous U (1)a 
scenario, A can be estimated as 

{A§)xx~>l fi - f 'K\f<-fi\)<DA)> (6-4) 

since the mass difference is given by (fa — fj) (Da), where fa is the anomalous U(\)a 
charge of F{. Here, the reason for the appearance of the coefficient W is that 



*^ The large SUSY breajking scale can make it possible to avoid the flavour changing neutral 
cur rent (FCNC) problem, H'Elf but in our scenario this is not the case, because the anomalous 
U(1)a charge of the Higgs H is inevitably negative, which forbids the Higgs mass term at the tree 
level. 

**' The CP violation parameter €k gives constraints on the imaginary part of (<5i2)xy that are 
approximately one order more severe than those it places on the real part. Here we concentrate 
only on the constraints from the real part of K°K° mixing, since under the other experimental 
constraints on the CP phase originating from SUSY breaking sector, which are mainly given by 
the electric dipole moment, we may expect that the CP phases are small enough to satisfy the 
constraints from the imaginary part of the K°K° mixing. 
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the unitary diagonalizing matrices are given by 

( 1 \\h-f 3 \ 



(6-5) 



In our scenario, the anomalous U(1)a charge of 5i is the same as that of 52; i.e., 
the sfermion masses of 5% and 52 are almost equal. This weakens the constraints 
from these FCNC processes. This is because the constraints from K — K° mixing 
and CP violation on the product {5\2)ll x {$12) rr are much stronger than those 
on {5i2) 2 LL or (5i2) RR , as shown in Eqs. (6T) and (6-2). Therefore, suppression 
of {Ai^)rr makes the constraints much weaker. Because the constraints from the 
mixing (and the CP violation) become weaker, as discussed above, we have 
a larger region in the paramter space where lepton flavour violating processes like 
[i — > &y are appreciable. Actually, if the ratio of the VEV of Da to the gaugino mass 
squared at the GUT scale is given by 

(Da) 



'1/2 

then the scalar fermion mass square at low energy scales is estimated as 

m% ~ fiRM? /2 + r) F M? /2 , (6-7) 

where rjp is a renormalization group factor. Therefore, in our scenario, Eq. (6-2) for 
(5 12 )ll becomes 



~ A ^~!f Ml '\ = A to-£> R (6-8) 

2 / (VD L + ^4^ J R) 1/2 M 1/2 (GeV) \ , 
< 4.0 x 1(T 2 I — 2 ; tlfl L j , (6-9) 



which can be rewritten 

M 1/2 > 1.25 x 10 4 A-^-^— -(GeV). (6-10) 



Though the main contribution to LS^RR vanishes, through the mixing in Eqs. (4-4) 
and (4-5), (Sfy) rr is estimated asfcil 

(^) M ~Ai A '^-f l ' B , (641) 
vd r +^\R 

where the mixing A 2 is different from the naively expected value 1 = A^ 1- ^ 1 . From 



M 1/2 >1.8xlO ^-^M (6 , 2) 



Eq. (6T) for J (5i 2 ) l l(S 12) RR, the constraint on the gaugino mass Mx/ 2 is given by 

, A 1 - 75 fi(^i - V > 2 ) 



*' We thank S. Yamashita for pointing out the contribution from the normalization factor of 
the main mode. 
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Mx/2 > 1.6 x 10 



3W1 



1/2 



VE R + 



-(GeV). 



(6-13) 



Taking the reasonable values f/'i = 5, ip 2 = 4, ?xd l ~ ry^^ ~ 6 and r]E R 
lower limits of the gaugino mass are roughly given as in Table I. 



0.15, the 



Table I. Lower bound of gaugino mass M\i 2 at the GUT scale (in GeV). 



7? 


0.1 


0.3 


0.5 


1 


2 




15 


38 


53 


73 


86 




60 


150 


210 


280 


350 


(5f 2 )Kil 


370 


260 


210 


150 


110 



Note that when R_< 0.5, the /i — > e^f process gives the severest constraint intheae 
FCNC processes. ElP We conclude that the lepton flavour violating processes!^ 'EHP 
might be seen in the near future. 

The reason for the suppression of (A® 2 )rr is that the anomalous U(1)a charge 
of 52 becomes the same as that of 5i, because the fields 5i and 5 2 originate from 
a single field, This is a non-trivial situation. The massless mode of the second 
generation 52 = $i(10, 5) + A 5 / 2 <p3(16, 5) has Yukawa couplings through the second 
term A s ' 2 !p3(16, 5). However, for the SUSY breaking term, which is proportional to 
the anomalous U(1)a charge, the contribution from the first term dominates that 
from the second term. This results in degenerate SUSY breaking terms between the 
first and the second generation. It is obvious that the twisting mechanism in E$ 
unification plays an essential role in realizing this nomjtrivial structure. Note that 
such a structure is realized only when (§1, 82) = (1, l')pl in which bi-large neutrino 
mixing angles are also realized. It is suggestive that the requirement to reproduce 
the bi-large mixing angles in the neutrino sector leads to this non-trivial structure, 
which suppresses the FCNC processes. till In this way, such a non-trivial structure is 
automatically obtained in the Eq model, which is much different from the situation 
for the SO(10) model, in which the condition can be satisfied only by hand. 



*' The case (1, l',2') cannot realistically yield the large mixing angle indicated by atmospheric 
neutrino experiments. 

**' We should comment on the D-terpi contribution to the scalar fermion masses. Generically, 
such a D-term has non- vanishing VEVE3 when the rank of the gauge group is reduced by the 
symmetry breaking and SUSY breaking terms are non-universal. In our scenario, when the E§ 
gauge group is broken to the S'O(IO) gauge group, the D-term contribution gives different values 
to the sfermion masses of 16 and 10 of S'O(IO). This destroys the natural suppression of FCNC in 
the E% unification. However, if SUSY breaking parameters become universal for some reason, the 
VEV of D can become negligible. Actually, the condition m| = m| causes the VEV of the D to 
be greatly suppressed. Therefore, in principle, we can control the D-term contribution, though it is 
dependent on the SUSY breaking mechanism. 
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§7. Discussions and summary 

In this paper, we examined an Eq unified model in which bi-large neutrino mix- 
ing angles are realized. A noteworthy fact of such GUT model with the anomalous 
U(1)a framework, is that once we fix the charges of all the fields of the model, all the 
hierarchical scales, the symmetry breaking scales at high energy and also the hier- 
archical structure of the Yukawa couplings, are determined without any ambiguity. 
The only exceptions are the SUSY breaking scale and the electroweak breaking scale, 
which we here adjusted from the experimental W masses. Even if the SUSY breaking 
scale is introduced by hand, we have to explain why the SUSY Higgs mass param- 
eter fi is around the SUSY breaking scale (the \i problem). One possible solution 



of the // problem has recently been examined in Ref. 34). Here we summarize the 
essence of the finding there. The SUSY Higgs mass, which is forbidden by the SUSY 
zero mechanism, can be induced when SUSY is broken. Thus the parameter \x must 
be proportional to a SUSY breaking parameter, and its coefficient is determined by 
anomalous U(1)a charges. Let us introduce the GUT gauge singlets S, with positive 
s, and Z, with negative z, with the mass term \ S+Z SZ in the superpotential. Since S 
has positive charge, it has vanishing VEV in SUSY vacua (see the Appendix). When 
SUSY is broken, generically a tadpole term X s AS (A is a SUSY breaking parameter) 
is induced in the SUSY breaking potential Vsb- As a result, the S field develops a 
non- vanishing VEV as 

(S) = \- s - 2z A. (7-1) 

Using this VEV shift, we generically find the \i term to be proportional to the SUSY 
breaking paramter A. In our Eq scenario, introducing the superpotential 

W = \ s+rp+2h S$H 2 , (7-2) 

the SUSY Higgs mass \x is obtained as 

fi ~ A 2(ft-*)+A(^) A (7 . 3) 

Therefore, if 

-l<2(fc-s) + i(0-0)<l, (7-4) 

the \i parameter becomes naturally around the SUSY breaking scale. Moreover, 
the Eq gauge singlet fields S and Z can be identified with composite operators, for 
example, we can take Z ~ CO or Z ~ <P<P. 

In our Eq case, the minimal field contents are, in addition to (9, three matter 
multiplets, #i(27), the pair of Higgs fields ^(27)^(27) and the pair of Higgs fields 
C(27) and (7(27), which are needed for the breaking Eq -> 50(10) -> SU(5) -> the 
standard gauge groups, together with an adjoint field A(78), which also provides a 
natural D-T splitting mechanism, as explained in the context of the SO (10) model 
in separate papers, t and leaves light Higgs doublets H. Among those minimal 
contents of matter and Higgs fields, we have nine charges, (i(ji,ip2,tp3), (<j>, 4>), (c, c), a 
and h, which determine the main features of the mass matrices of quarks and leptons. 
First, the CKM mixing angle almost fixes the charges of the matter fields ipi = 
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(n + 3, n + 2, n), and the doublet Higgs /i = — 2n, and in order to get bi-large mixing 
angles for the neutrino sector, we need a constraint on the charges of the Higgs fields, 
r = 1/2, i.e., c — 5=0 — + 1, and also we have the constraint <p — <j) = 2n — 10 — I 
(—2 < I < 2) in order to give the proper neutrino masses. If we choose the charge as 
a = —2 in order to insure that the GUT relation between the masses of down-type 
quarks and charged leptons acts only for third generation, there remain only three 
degrees of freedom. Moreover, it may be possible to build DT splitting models in 
which the light Higgs can be identified with the components of <P or C. Actually, 
it is naturally realized that <P can play the role of H, so in that case, <fi = h.c3 
This implies that there are now two degrees of freedom. If we further impose the 
condition of solving the n problem ( |7-4| ), there remains just one degree of freedom. 
It is quite suggestive that there is a set of charge assignments that can satisfy all the 
above conditions. Actually, if we take n = 2, (p = h = —4, <p = 0, c = —4, c = —1, all 
the above conditions are satisfied for I = —2. The charge assignment n = 2, 4> = h = 
—4, cf> = 3, c = —6, c = 0, 1 = 1 is quite interesting, because the composite operator 

can even play the same role as the field when £ 2 ~ A; that is, (p<P) ~ A = £ 2 . 
In this case, we have the minimum model in which there are *Zq, ^3) C, C 

and A, where all the charges are uniquely determined. 

What is interesting in the Eq unified model is that the condition for suppression 
of FCNC is automatically satisfied. The essential point is that the first and second 
generation fields of 5 have the same anomalous U{\)a charge because these fields 
originate from a single field <Zq. 

The aspect of the family structure that has recently been made clear by the 
neutrino experiments gives a guide to investigate the origin of the family. The 
scenario discussed here is quite impressive, and it leads us to expect that we may 
find "the real GUT" in the near future. 
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Appendix 



In this appendix, we explain how the vacua of the Higgs fields are determined 
by the anomalous U(1)a quantum numbers. 

First, we show that none of the fields with positive anomalous U(1)a charge 
acquire nonzero VEV if the Froggatt-Nielsen (FN) mechanism acts effectively in the 
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vacuum. Let the gauge singlet fields be (i = 1, 2, ■ • • n±) with charges zf , when 
z~l > and < 0. From the F flatness conditions of the superpotential, we get 
n = n+ + n_ equations plus one D-flatness condition, 

^ = 0, D A = g A [y2z t \Z i \ 2 + f\ =0, (A-1) 



{^Zi\Zi 



where £ 2 = g ° 9 ^ A (= X 2 Mp). At first glance, these look to be over-determined. 
However, the F flatness conditions are not independent, because the gauge invariance 
of the superpotential W leads to the relationtH 

SW 

i^rZiZi = 0. (A-2) 

Therefore, generically a SUSY vacuum with (Zj) ~ Mp exists (Vacuum a), because 
the coefficients of the above conditions are generally of order 1. However, if n + < n_, 
we can choose another vacuum (Vacuum b) with i^t} = ^> which automatically 

satisfies the F-flatness conditions = 0. Then the (Z^j are determined by the 

F-flatness conditions = with the constraint ( |A-2j ) and the L>-flatness condition 

i 

D A = 0. Note that if A < 1 (i.e., £ < 1), the VEVs of Zf are less than the Planck 
scale. This can lead to the Froggatt-Nielsen mechanism. If we fix the normalization 
of U{1)a gauge symmetry so that the largest value in the negative charges z~[ 
equals -1, then the VEV of the field Zf is determined from Da = as ~ 

which breaks U{\)a gauge symmetry. (The field Z-j - was introduced in the previous 
section as 0.) Other VEVs are determined by the .F-flatness conditions of Zf as 
ZT > ~ A _z i , which is shown below. Since (ZJM = 0, it is sufficient to examine 



the terms linear in ZA in the superpotential in order to determine \ Z i j. Therefore, 
in general the superpotential can be written 

W = Y,W Z+ , (A-3) 



\ 3 j,h J 

Ti-\- / n_ n_ \ 

= E# E^T+E^ + - . ( A ' 5 ) 

* \ 3 j,k ) 

where Z = A 2 *Zj. The .F-flatness conditions of the Zf fields require 

X< (l + E^~ + "--j = °> (A-6) 



*' We thank H. Nakano for pointing out this relation. 
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which generally lead to solutions Zj ~ O(l) if these .F-flatness conditions determine 
the VEVs. Thus the F-flatness condition requires 

(Z j )~0(X-'i). (A-7) 

Here we have examined the VEVs of singlets fields, but generally the gauge invariant 
operator O with negative charge o has non- vanishing VEV (O) ~ A~° if the F-flatness 
conditions determine the VEV. 

If Vacuum a is selected, the anomalous U(1)a gauge symmetry is broken at 
the Planck scale, and the FN mechanism does not act. Therefore, we cannot know 
the existence of the U(1)a gauge symmetry from the low energy physics. On the 
other hand, if Vacuum b is selected, the FN mechanism acts effectively, and we can 
understand the signature of the £7(1) a gauge symmetry from the low energy physics. 
Therefore, it is natural to assume that Vacuum b is selected in our scenario, in 
which the U{\)a gauge symmetry plays an important role for the FN mechanism. 
This amounts to assuming that the VEVs of the fields Zf vanish, which guarantees 
that the SUSY zero mechanism acts effectively. 
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